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Area  of  frame 
IntegraUon  constants 
Bay  matrix;  see  Equation  [8] 

Faying  width  of  frame 
Eh^ 

12(1 -k-*) 

Young’s  modulus 
Frame  matrix 
Deflection  functions 
Deflection  functions 

(Ajp  +  AAV^Vi  +  o*/3^  sinh  Ajl  +  Vl  ~a^$^  sin  Ajfj  + 
2A  - 

“cTrl  (cosh  K  -  cos  K^l) 

Shell  thickness 

Moment  of  inertia  of  frame  section  about  radial  axis 

o  Ki  +© 

Center-tO'Center  distance  of  frames 
Clear  distance  between  frames 
LongitOdinal  bending  moment  in  shell 
Pressure 

Collapse  pressure 

Radius  to  middle  surface  of  shell 

Radius  to  center  of  gravity  of  frame  section 

Bay  matrix;  see  Equation  [11] 

2Aa  j 

Longitudinal  shear  stress  in  shell 
Distance  along  shell  generator 
Radial  deflection  of  shell 
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Differential  shear  function 

Poisson's  ratio 
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Functions  defined  by  Equations  [41],  [44],  [46],  and  [48l, 
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ABSTRACT 


Recent  developments  in  structural  research  on  circular  ring-stiffened 
cylinders  subjected  to  hydrostatic  pressure  have  indicated  that  refinements  of 
the  standard  strength  analysis  are  required  to  account  for  the  effects  of  discon¬ 
tinuities  such  as  a  heavy  frame  or  bulkhead  or  variations  in  shell  thickness.  A 
procedure  has  been  developed  for  computating  axisymmetric  stresses  near  these 
discontinuities  in  a  cylinder  in  which  the  thickness  of  the  material  does  not  vary 
between  a  pair  of  stiffeners  but  may  change  from  one  side  of  a  stiffener  to  the 
other.  The  results  of  this  analysis  indicate  that  stresses  higher  than  those  pre¬ 
dicted  by  a  standard  computation  will  usually  exist  near  these  discontinuities. 

A  method  of  reducing  this  effect  by  modifying  the  geometry  near  these  points  is 
also  presented. 


INTRODUCTION 

The  strength  of  thin-walled,  stiffened,  cylindrical  shells  against  axisymmetric  yielding 
is  usually  computed  from  an  elastic  analysis  of  Von  Sanden  and  Gunther.*  This  analysis  gives 
the  stress  distribution  in  the  shell,  and  colliqjse  is  assumed  to  occur  when  the  stresses  at  the 
middle  of  a  bay  exceed  the  .yield  point  of  the  material.  A  more  recent  analysis  developed  by 
Salerno  and  Pulos^  accounts  for  the  longitudinal  moments  resulting  from  the  end  pressure  to¬ 
gether  with  the  radial  deflection  of  the  shell  and,  therefore,  gives  more  accurate  stresses. 
Although  neither  of  these  analyses  indicates  any  difference  in  strength  between  bays,  it  has 
been  observed  experimentally  that  failure  by  axisymmetric  yielding  almost  invariably  occurs 
in  the  full-length  bay  nearest  a  bulkhead  or  heavy  frame.  To  determine  the  effect  of  end  con¬ 
ditions  on  the  stresses  in  the  shell,  a  modified  Salerno  and  Pulos  analysis  is  derived  which 
takes  into  account  variable  frame  spacing  and  size  and  changes  in  shell  thickness;  it  can  be 
easily  extended  to  variations  in  radius  between  adjacent  bays. 

To  offset  the  weakening  effect  of  the  heavy  member,  it  has  become  customary  to  reduce 
the  spacing  of  the  first  frame  from  this  member.  This  practice,  however,  merely  forces  the 
weakening  effect  of.  the  heavy  member  into  the  next  bay  without  appreciably  increasing  the 
collapse  pressure.  Since  it  is  usually  undesirable  to  reduce  the  length  of  more  than  one  bay, 
another  method  to  increase  the  strength  of  Uie  end  bays  is  proposed  in  this  report. 

ELASTIC  STRESSES  IN  FINITE  CYLINDER 

The  theory  of  Salerno  and  Pulos  assumes  an  infinite  closed  cylinder  composed  of  a 
linearly  elastic  material  of  uniform  thickness  and  radius  reinforced  by  equally  spaced,  equal¬ 
sized  frames  and  subjected  to  hydrostatic  pressure.  A  typical  bay  of  this  cylinder  illustrating 

*R«(«r«ic«a  are  liatad  on  pmgo  24. 
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- -  - = - -  ^  Cylindtr  —  - -  ♦ - 

Figure  1  -  Typical  Bay  of  the  Salerno  and  Pulos  Cylinders 

the  initial  and  deformed  shapes  is  shown  in  Figure  1. 

The  deformation  defined  by  the  radial  deflection  y  as  a  function  of  z  must  satisfy  the 
differential  equation 


d*y  ^  12(1-./^  PR  d^y 
dz*  ■  2Eh^  dz^ 


{l-u/2) 


]• 


0 


[11 


where  E  is  Young’s  modulus,  u  is  Poisson’s  ratio,  and  the  other  terms  are  defined  in  Figure  1. 
Setting 


/3^ 


PR^ 

IeK' 


and 


PR^ 

IE 


this  simplifies  to* 

y*''  +  4a^/3^  y"+  4o^(y-y,)  *  0 


[21 


Note  that  y^  is  the  deflection  which  would  occur  if  the  stiffening  frames  were  not  present. 
One  form  of  the  general  solution  of  Equation^  [1]  or  [21  is 

y  =  i4j  cosh  K  ^z  Cos  +  ^2  ^1*  ^2®  ^3  ^  1*  ^2* 

+  sinh  K  ^z  cos  K^z  +  y^  [Sal 

where  *aVi-o*/3^  and  -a^l  +  The  derivatives  of  y  are  readily  obtained  in 

the  following  form 

y'«  (•^1^1  +  ^2^2!  ^1®  ^2*  [^2^1  “  ^1^2^  ^i*  ^2* 

+  (^3  Kj  +  ^4^1)  cosh^fjXcos  ^2*  +  (^3^1  “  ^4^2!  ^1®  **"  ^2* 

*Th«  MlttUon  which  follow*  i*  takaa  fioin  R*farMc*3. 
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y"  •  {A^(K\  -  K\t  +  jl  cosh  K\x  cos  /',»  +  -  2A^K^K^ 

sinh  sin  /fj*  +  [^<3(^1  -  K\)-  2A^K^Kj^  cosh  K  j»  sin  K^a 
+  [/I4  (f(\  -  +  2i43  K^K^]  sinh  AT,*  cos  K^t 

y'"m  {A^K^{K\  -  3^1 )  +  A^^ZK\  -  K^]  sinh  K^a  cos  Afj*  +  («'|  -  3A'f) 

+  -  3Aj  )\  cosh  /f  ,*  sin  j*  ^  IAJ(^ZK\  -  /T*)  +  A^K,  (/f  J  -  .W|)l 

oosh  K^x  coa  K^m  +  »  -  3A'^  +  A^K^kI  -  S/f?)!  sinh  K sin  K^x 

Let  _ 

11  -  (Ap  +  bk)  ( VTTo^  sinh  K^U  ainKJ)^-^  V^l  - 

(oosh  -  cos 

where  Ap  is  the  area  of  the  frame  cross  section*  and  6  is  the  faying  width  of  the  frame. 
Applying  conditions  of  continuity  and  equilibrum,  Equations  I3l  for  the  Salerno  and  Pulos 
cylinders  may  be  reduced  to: 


[3c] 


[3d] 


[41 


where 

/j(xl  -Vl  -  [oosh  Aj*  sin  ®osh  Aj(l-x)  sin  Ajx]  +V 1 

[sinh  KyX  cos  K^l-x)  +  sinh  Aj  (f-xl  cos  Ajx] 

2  [sinh  Aj*  sin  K^t-x)  ~  sinh  K^l-x)  sin  Aj^d 


[6al 

[6b] 


-  }f\  -  a [cosh  Aj*  sin  K^l^x)  +  cosh  A  sin  A -Vl  + 

2a* 

[sinh  Aj*  cos  A/t-x)  +  sinh  Aj(f-  x)  cos  Ajx] 

•OraatM  •ccuracy  i*  obUUi«<l  by  Ming  lor  Ap  tho  voluo  obUinod  by  mulUplylng  th.  tru«  of  tho  ooctton  by 
R/Rp  whoro  Rp  l«  the  rodioo  of  tho  contor  of  groTity  of  tho  fromo. 
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/»00 

1  (  r 

1a  ■ -  -  -  2  4  [sinh  K.x  sin  K  Jl-st)  -  sinh  K^{1-  ai)  sin  K +  Vl -m*P* 

2rt* 


[oosh  K^(l-x)  oos  -  cosh  K^x  cos 


[6d] 


In  an  actual  structure,  the  presence  of  bulkheads  or  heavy  frames  precludes  the  assump¬ 
tion  of  one  deformation  pattern  for  all  bays  as  assumed  in  the  Salerno  and  Pulos  analysis. 
However,  an  approximate  relationship  can  be  found  by  computing  y  and  its  derivatives  at  a 
distance  from  the  bulkhead  in  terms  of  the  conditions  at  the  bulkhead  and  by  assuming  that 
these  will  approach  the  values  computed  by  a  Salerno  and  Pulos  analysis  as  the  distance  in¬ 
creases.  Basically,  symmetric  deflection  assumptions  of  Salerno  and  Pulos  are  replaced  by  the 
matrix  equation 


l^lo.SAP 


[7] 


where 
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y' 
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II 

2«* 


2o» 
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[81 


The  subscript  0,  i  denotes  that  |61  is  to  be  evaluated  at  the  bulkhead  or  at  a  frame (x^O)  in 
the  tth  bay  from  the  bulkhead,  and  the  subscript  0,  S&P  denotes  values  at  x  -  0  in  the  infinite 
cylinder.  Since  Equation  [7J  is  rarely  exactly  satisfied,  i8lo,(  be  taken  such  that  it  is 
most  nearly  satisfied  by  a  “least-squares*’  approximation  consistent  with  the  bulkhead  con¬ 
ditions. 

The  basic  problem,  then,  is  to  determine  in  terms  of  the  bulkhead  condition  IBIq  j 
so  that  Equation  [7]  and  a  least-square  estimate  can  be  used  to  determine  If 

Equations  [3]  are  evaluated  at  x  =  0,  they  can  be  reduced  to 

-  y(0)  - 


*W«  could  ■■  ••tily  have  used  • 
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[9a] 


[9bl 

'  2. 

'  'Vi  ' 

(y(0)  -  y.) 

v'l  - 

1 

fy'(0)\l  1+2cPB^ 

1/  ' 

[9c] 

“F 

'  '  2-»  Y 

1 

Yy'(0)W  l-2a^B^ 

V  1 

[9d] 

°  2 

1  \  ic? 

Ayi-a^^^  1} 

Now,  substitution  of  Equations  [«!  into  Equations  [3],  which  will  yield  conditions  at  any  point 
of  the  bay,  gives 


^  '^•’en  a.  <  f. 

where  [51^^  .  is  defined  by  Equation  [11]  on  page  6. 

The  identity  equation  a  y^,  which  holds  throughout  the  bay,  is  included  in  Equation  [10] 
for  reasons  of  conformableness  which  will  be  apparent  later.  Replacing  >  by  1,  and  letting 
1  a  1  in  Equation  [10]  gives 

^  [121 

now  denotes  the  boundary  conditions  at  the  end  of  the  bay  away  from  the  bulkhead.  It 
is  possible  to  assume  that  |B|^  ,  -  l^l/.  saP  and  to  solve  immediately  for  IBI^  ,.  This  will 
not  give  a  very  good  approximation  because  f ,  is  usually  too  small.  In  order  to  obtain  |B| 
at  a  peater  distance  from  the  bulkhead,  |e|  must  be  found  on  the  other  side  of  the  frame.  *The 
conditions  on  both  sides  of  a  frame  are  shown  in  Figure  2. 
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Assuming  that  the  faying  width  b  remains  straight,  four  equations  can  be  obtained  from 
considerations  of  continuity  and  equilibrium.  Obviously 


Equilibrium  of  vertical  forces  yields 
V^d^-V.^.m-Pb^E 


[13al 

[13b] 


■  Vidi^  +  y/  + 1  (OV 

/Ap  +  bhy 

IJ 

2R 

1  R  i 

1" 

1  2  )\RI 

Moment  equilibrium  yields  ^3^ 

jr  1  H’f'Ti') 

»,  (I, )-«,*,«» -2  •'it .(»'] -S'/ ly — —  ■“ 


where  Ip  *  is  the  moment  of  inertia  of  the  frame  about  the  vertical  (radial)  axis.  It  should  be 
noted  that  h  may  be  selected  from  either  bay  as  long  as  bh  represents  the  total  area  of 
frame  and  faying  width  of  shell  and  Ip  +  (6^4/12)  represents  the  total  moment  of  inertia  about 
the  radial  centerline.  The  moments  and  shears  are  obtained  in  terms  of  y"  and  y*"  by  the 
relations 


-  Dy'r 

[14] 

-  Dy'" 

[15] 

where  D  « - for  a  cylinder. 

12(1  - 

Substituting  these  expressions  into  Equations  [13]  and  rearranging  gives; 

-  leiu 


where 


R 

*Graat«r  accuracy  is  obtained  by  using  for  Ip  the  value  obtained  by  eniltiplying  the  true  /  of  the  section  , 


T 


t' 
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The  double  subscnpt  1 4-  1  indicates  a  frame  between  i  and  » 4-  1  bays  from  the  bulkhead. 
Equation  [161  includes  the  identity  equation 


J  4-  1 


5 


(y.V 


and  allows  the  shell  parameters  (except  radius)  to  change. 
Substitutiori  of  Equation  [111  into  Equation  [161  gives 


1+1 


l/.i 


.) 


[181 


By  repeated  application  of  Equation  [181  we  may  obtain 


+  fn 


|S|,, 


i+m— 1 


X  . 


\r 


1.1+ 1 


X  X  \B\ 


O.i 


[191 


For  t  B  1  Equation  [  19]  may  be  shortened  to 


-If’SU*.,.  "Ifilo.i  1201 

where  IF,  j  is  the  product  of  the  [Fj  and  |S1  matrices  in  the  order  of  Equation  [19]  for 

all  frames  and  bays  between  the  (m  +  1)  bay  and  the  bulkhead.  Now  if  the  product  for  matrix 
IF,  S|^^j  1  with  at  least  two  consecutive  identical  bays,  Equation  [7l  will  be  very  nearly 
satisfied  by  i  «  m  +  1.  Almost  all  problems  of  practical  interest  can  be  carried  far  enough  to 
satisfy  this  condition,  and  it  will  be  assumed  hereafter  that  |F,  i  does  in  fact  terminate 

at  such  a  point.  Substituting  Equation  [20]  into  Equation  [71  gives 


01+  1 , 1 


IBIo.i 


[211 


As  previously  noted  this  is  not  exactly  true,  but  a  least-squares  process  can  be  used  to 
estimate  two  of  the  bulkhead  conditions  if  the  other  two  are  given.  Any  pair  may  be  used  as 
the  knowns,  buty  and  y"  will  be  used  here  because  they  are  the  most  readily  evaluated  at  a 
bulkhead  at  pressures  sufficiently  high  to  cause  yielding  at  the  bulkhead.**  Setting 


1 


*n  now  indlcatM  Uie  number  of  framea  from  any  atartlng  bay  i. 

**See  Appendix  A. 
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ud  B«gl«otiBg  the  identity  equation,  Equetioo  [21]  redueee  to 

«1«  *  PfV  ^  m  ^)  [22al 

Ojtt +/32  »  +  Vj  -  0  t22bl 


"3“  +  »  +  ya  -  0 

®4“  +  ^4*+y4-0 

The  method;  of  least^squares  shows  that  u  and  v  are  given  by 

lapiPy-ttaylft^ 


[22c] 

(22dl 

[23] 


.  -  [841 

-  (iapy 

Combining  u  and  v  ■with  tW.known  values  y,  y'72o*  ,  and  gives  the  complete  column  matrix 
|8|q  j ,  and  the  column  matrix  |6|q  ^  at  any  other  frame  is  then  found  by  multiplication  by  the 
IF,  S|^.  j  matrix.  Equation  [10]  then  completely  determines  the  deformations,  and  hence  the 
stresses,  in  any  bay. 


OPTIMUM  DESIGN  FOR  END  BAY 

The  analysis  of  the  previous  section  is  applicable  to  ui  axially  symmetric  structure 
composed  of  circular  cylindrical  bays  with  a  common  axis  and  radius.  The  length  of  bay, 
thickness  of  shell,  and  Young’s  modulus  may  vary  from  bay  to  bay;  the  size  of  stiffeners  may 
also  change. 

In  actual  practice,  however,  all  these  parameters  are  usually  constant  except  for  iso* 
lated  interruptions  such  as  bulkheads,  heavy  frames,  and  conical  reducer  sections.  Analyses 
made  on  several  cylinders  in  Uie  range  of  current  submarine  design  indicate  that  there  is 
usually  a  bay  near  any  bulkhead  or  other  variation  of  geometry  which  has  larger  elasUc  circum* 
ferentiai  membrane  stresses  than  those  bays  farther  from  the  bulkhead.  The  following  design 
method  has  been  developed  to  eliminate  this  weakness. 

For  this  design  procedure  a  long,  ring*stiffened  cylinder  of  uniform  geometry  composed 
of  a  ductile  material  and  terminating  at  a  relatively  rigid  bulkhead  or  other  discontinuity  is 
assumed.  It  is  desired  that  the  stress  pattern  caused  by  the  bulkhead  shall  not  cause  any  bay 
to  be  weaker  than  a  bay  more  distant  from  the  bulkhead.  It  is  apparent  that,  if  Equation  [21] 
is  exactly  satisfied  for  some  m,  then  all  subsequent  bays  will  have  the  Salerno  and  Pulos^ 


10 


stress  distribution.  To  sstisfy  Equation  two  additional  parameters’ are  required.  If 
they  are  selected  from  the  geometry  of  the  first  bay  and  first  frame,  and  Equation  [21]  can  be 
satisfied  for  m  -  1,  the  stress  pattern  of  the  second  and  all  subsequent  bays  will  be  identical. 
The  length  of  the  first  bay  and  the  area  Ap  of  the  first  frame  are  selected  as  the  parameters 
to  be  determined.  Since  stresses  are  nonlinear  function^  of  pressure,  this  determination  should 
be  made  for  P  ■  P^,  the  expected  collapse  pressure. 

It  is  more  convenient  now  to  select  the  origin,  e  ■>  0,  of  the  first  bay  at  the  first  frame, 
while  the  origin  for  the  second  or  typical  bay  will  remain  as  it  was.  The  nomenclature  for  this 
case  is  shown  in  Figure  3. 


Figure  3  -  Nomenclature  for  End  Bay  Design 
At  the  right  of  the  frame  we  set 


From  Equation  [16] 


iSl  o  2  ”  |Plo,S4/> 


[25] 


[26] 


where  ^is  identical  to  |  with  the  second  and  fourth  columns  multiplied  by  (-1). 

Since  y'g^p  (0)  »  0,  substituting  Equation  [25]  into  Equation  [26]  and  rearranging  yields 


Vi(0)  -  ysfc/»(0)  =  0 


{27a] 


y,'(0) 


— = — P  0 


[2Tb] 


11 


[27cJ . 


yr(0) 

~2^ 


— - “IT"  ^ 

«v 


yr^o)  ys;v(o>  +  w 

- ; - 2  (6«i  — —  (0)  -  2yslp  (0)  ^  2  {ha)y^ 

2rf  2fl^ 

It  is  now  convenient  to  define  5  and  AV  as  follows: 


5  ■»  "'  —  a  ■'  ''  '  ■■"  '  ■■■  I  ■ 

it^D  it?D 

whence  Equations  [27cl  and  [27d]  may  be  written  as 


yr(0>  ysLp(^) 

- m.(bd)B 

2«^  2a* 


yrm  ysApfo) 

- -  =  -  2  a 

2a*  2a* 

The  difference  between  Equation  [10]  for  Bay  1  and  Equation  [10]  for  Bay  2  is 


[27dl 


[28] 


[27c'] 

[27d'] 


*,2 


[20] 


since  under  the  assumptions  of  this  section  is  the  same  for  all  bays. 


[27] 


l^lo,  1  "  1^1' 0,2 


0 

0 

(6a)S 

25 

0 


From  Equations 


[801 


and,  for  ar  -  fj,  Equation  [29]  becomes 


0 

0 

(&o)5 

25 

0 


[311 


Now,  since  |S|^  ■  and  |61^  2  completely  known  from  the  geometry  and  Equations  [11]  and 

[5],  respectively,  and  since  there  must  exist  two  boundary  relations  on  Equation 

[31]  represents  four  equations  in  the  two  unknown  boundary  conditions  and  the  two  parameters 
f  j  and  S.  Equation  [28]  defines  5  as  a  function  of  the  difference  between  the  shear  on  one  side 
of  the  first  frame  and  that  on  the  other  side.  This  indicates  that  the  area  of  the  first  frame 


•This  subject  is  discussed  in  Appendix  A. 
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must  be  altered  to  maintain  the  assumed  oonditionb.  Since  the  S&P  frame  is  in  equilibrium 
except  for  the  load  \V,  a  corresponding  A  is  computed  from  the  equilibrium  expression 


and  Equations  [Sal  and  [281  as 


®y2(0) 


[321 


assuminii;  a  solution  to  Equation  [31]  exists  then  Equation  [8  ll  and  the  two  known  boundary 
relations  will  completely  determine  the  geometry  necessary  for  the^ondition  that  the  second 
and  all  subsequent  bays  will  have  identical  stress  patterns  at  the  collapse  pressure.  This  is 
developed  in  detail  in  Appendix  A.  A  sample  solution  is  presented  in  Appendix  B. 

To  determine  the  strength  of  the  first  bay,  Equation  [29]  is  written,  after  substitution 
of  Equation  [301,  as  follows; 


yi(*)  - 

[34al 

yi{x)-y{{x) 

=  -  5  ^2  (3!) 

[34bl 

[34c] 

2a^ 

y'"{x)-y^"{x) 

e  -  5 g^{x) 

[34d] 

2a^ 


where 


cosh  K,  *  sin  K,  x  sinh  K,  *  cos  x  sinh  K^x  sin  K^x 

*  — - - — 


[35  a] 


*Th«  MtMd  •ddiUonal  requirad  will,  of  oouroo,  be  given  by 
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9{(ah .  2  8inh  /if}(» 
^ 

sioh  K^a  ooB 

^  / 


I /if}  (t  ain  IfjS  /ooah  ^f} » 

: . .  . . +  (>^8  ’  '  - -  - 

1-«V  \  ¥l+a^$* 

'aA 


oodilfiaain/fj'  siiih  sooa  Af,  «  ^ 

(*)  ■  — — >'■!■— - •’—■  + - : -  ^ — —  +  bti  (coBh 

VT+«^  ' 

flf*/8*  \ 

-y—  •  sinh  K.  a  sin K^a) 

Vl-u*fi*  ‘  / 


(35b] 

«C08  ^2  * 

[35c] 


■  —  .  2  ooah  Ki  aooa  -r  -■•■  -  sinh  Jf.  a  ain  a 

2o* 

'1-W02  l4.2a%2 

sinh  Af}  a>  coa  Af»  a; -  — cosh, 


.  K/izffV 


i-  *_  _____ 


Af}  a  sin  ^ 


[35d] 


For  cylinders  of  interest  here,  it  will  be  true  that  <  yj  (/,>,  yj"<ij>> 

and  1}  < -~.  If  this  is  so,  then  y j  (*)  >  0  for  0<  *  ^  f,  and  y,  {a)  >  0  near  the  middle  of  the 

bay.  Thus  the  maximum  longitudinal  stress  near  the  middle  of  the  first  bay  will  be  somewhat 
higher  than  that  of  a  typical  bay.  However,  the  average  and  also  the  maximum  circumferential 
membrane  stresses  of  the  first  bay  are  leas  than  those  of  all  subsequent  bays.  Hence  the  first 
bay  should  withstand  a  higher  pressure  than  the  typical  bays. 


DISCUSSION 

The  analysis  of  this  report  is,  of  course,  limited  to  small  deflections  in  axially  sym¬ 
metric,  thin-walled  cylinders.  Therefore,  cylinders  with  the  ends  designed  by  the  method  out¬ 
lined  here  which  fail  by  an  axisymmetric  shell  yield  should  have  the  maximum  obtoinaMe 
strength.  It  is  logical,  also,  to  expect  such  cylinders  which  fail  by  shell  instobility  to  have 
at  least  as  high  a  collapse  pressure  as  can  be  obtained  with  any  other  end  design. 

In  applying  the  analysis  of  this  report,  there  are  a  few  pointo  to.be  remembered.  This 
analysis,  which  is  analogous  to  the  Salerno  and  PulOs  solution  for  the  infinite  stiffened  cylin¬ 
ders,  gives  results  which  are  nonlinear  functions  of  pressure.  Familiarity  with  the  Salerno  and 
Pulos  solution  should  not  be  allowed  to  lead  one  to  believe  that  this  nonlinearity  is  negligible, 
since  the  variation  in  deflection  near  a  rigid  or  nearty  rigid  bulkhead  may  easily  be  five  or  six 
times  that  in  a  bay  remote  from  the  bdkhead,  thus  greatly  magnifying  the  "beam-column"  effect. 
In  addition,  at  pressures  near  collapse  there  will  be  a  change  of  boundary  conditions  due  to 
higher  local  bending  stresses  near  the  bulkhead  which  produce  a  "plastic  hinge"  at  a  piessniu 
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w«ll  b«low  that  which  causes  collspss  of  the  oylisdsr.  Both  of  these  noBUnesr  effects  are 
even  more  important  where  an  outward  negative  deflection  is  expected,  such  ss  at  the  juncture 
of  a  cylinder  with  the  large  end  of  a  cone.  For  these  reasons  caution  must  also  be  exercised 
in  interpretation  of  experimental  strain  data  firom  a  cylinder  designed  by  the  methods  presented 
hers,  since  the  data  taken  at  low  pressures  will  always  show  greatsr  strains  near  the  middle 
of  the  second  bay  than  jn  any  other  bay.  It  is  only  whsn  the  pressure  nesre  the  collapse 
pressure*  that  the  strains  in  the  second  bay  will  approach  those  of  subsequent  bays. 
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AI^ENOIX  A 

BOUNDARY  CONDITIONS  AND  SPECIFIC  PROBLEMS 


Th«  U80fnlii«88  of  thi8  an8ly8i8  in  de8iKn  18,  of  course,  limited  by  the  ability  of  the 
deoigner  to  eetimate  boundary  oonditiono.  Although  it  is  not  necessary  that  the  boundery 
condition  be  known  absolutely,  an  attempt  should  be  made  to  limit  the  sum  of  the  absolute 

values  of  the  errors  in  each  of  the  pair  of  ^,(1 ; - i —  which  are  used 

*  *  «  2  01^  2  or 

to  5  percent  of  .  AlUiough  it  will  usually  be  necessary  for  the  designer  to  exercise  judg* 
ment  in  selecting  proper  boundary  values,  a  few  more  common  cases  will  be  discussed  here. 

In  all  the  examples  that  follow  it  will  be  assumed  that  the  material  is  ductile  and  of  the  ideally 
plastic  type,  In  other  words,  it  has  a  3tress*strain  curve  which  can  be  approximated  by  a 
linear  elastic  curve  terminating  in'a  well-defined  plateau-type  yield.*  It  will  bh  assumed  that, 
for  these  cases,.only  Ij  and  Af  ^  can  be  varied  to'attain  the  maximum  collapse  pressure. 

SYMMETRICAL  HEAVY  FRAME  OR  BULKHEAD 

By  a  bulkhead  is  meant. a  structure  for  which  it  may  be  assumed  that  all  deflections  are 
zero,  i.  e.,  Viili)  -  0.  It  is,  therefore,  only  a  special  case  of  a  heavy  frame  located  so  that 
it  may  be  assumed  to  be  symmetrically  loaded.  For  the  heavy  frame  there  are  two  possibilities 
to  consider;  (1)  a  “plastic  hinge"  is  formed  in  the  shell  at  the  frame,  or  (2)  partial  or  no  yield¬ 
ing  occurs  near  the  frame  prior  to  collapse.  Since  in  case  (1)  only  local  yielding  will  occur, 
it  is  assumed  for  the  First  case  that  the  tangent  to  a  generator  will  rotate  at  the  bulkhead  or 
heavy  frame  i.  e.,  y  j'  For  the  second  case  it  is  assumed  that  no  rotation  occurs  and 

y{<it)~o. 

For  the  first  case,  then,  two  conditions  at  the  frame  are  known.  The  first,  which  comes 
from  equilibrium  (see  Equation  [13o])  and  symmetry  (i.e.,  equal  deflections  and  shear  forces  on 
either  side  of  the  frame),  is  _  _ 


yr(ix) 


+  6^0 


Here  the  subscript  0  refers  to  the  heavy  frame.  The  second  condition  is  based  on  the  maximum 
strength  of  the  material  at  the  frame  which  occurs  under  the  assumed  equilibrium  stress  disr 
tribution  of  Figure  4.  It  is  —  P^R 


yt 


R  2  hr 


•Other  type*  of  material  are  beyond  the  acope  of  this  report 
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|— ^2 — ^ 


Figure  4  -  Assumed  Stress  Condition  at  Formation  of  Plastic  Hinge 


If  shear  stresses  are  neglected,  the  Hencky>Von  Mises  criterion  of  yielding  gives: 


<Tj- 


1 

2(1- v  +  v^) 

R 

1 

EViih) 

2(1->/+v2) 

R 

[38al 


[38bl 


where  Is  the  yield  stress  in  a  one>dimensional  stress  field  (assumed  to  be  the  same  in 
tension  and  compression).  Since  a  j  and  are  nonlinear  functions  of  it  is  easier  to 

make  a  preliminary  estimate  of  and  then  to  compute  these  stresses  directly.  Stresses 
so  obtained  will  usually  be  sufficiently  accurate. 

Combining  Equations  [86],  [34al,  [34bl,  [84d},  and  [5]  with  ar  »  Ij  gives: 


Combining  Equation  [371  with  Equations  [34c]  and  [5cl  gives: 


[40] 


where  T  ■  -2ho^/P^  R.  EliminaUon  of  5  from  Equ&tions  l89l  and  [40]  yields 

r  'a"-*  1  '*'• 

- T - 


[41] 


The  smallest  value  of  Ij  which  satisfies  Equation  [41]  is  the  desired  solution.  Since  for  f| 
near  /,  is  an  approximately  linearly  decreasing  function  of  1|,  Equation  [41]  can  usually  be 
solved  most  easily  as  follows: 

Step  1:  Compute  (1). 

Step  2:  If  0|(1)  <  0 ;  compute  (0.9  J). 

If  <^|(/)  >  0;  compute  (1.11]. 

Step  3:  Compute  by  a  linear  interpolation.* 

If  4>x  [0 


If  <f>^  (0  >  0,  then 


r  0.1  ^j(0  1 

1^‘  ^,(f)-0j(O.0oJ 
0.1  ^,(0  “I 


[42a] 


l>  I 


[42b] 


This  process  may,  of  course,  be  repeated  for  more  accuracy.  Xor  the  case  of  partial  or  no 
yielding  at  the  heavy  frame,  we  assume  y{ili)  -  0.  It  follows  from  Equation  [84b]  that 
Equation  [40]  can  be  replaced  by 

Af 


[43] 


*M6r«  accuracy  may  be  obtained  by  computing  an  Interinadlate  value  of  tft^  and  paaalng  a  circle  throu^  the 
three  valuea  graphically. 
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$ 


■  .'iwr, 


Thm  th*  •xprassioD  analogous  to  ^Equation  [41]  is 


y,«,)  ~1  it. 


[441 


The  solution  is  by  the  same  method  as  that  used  for  Equation  [41].  The  correct  value  of  f| 
is  then,  of  course,  the  smaller  of  Uie  two  values  saUsfying  Equations  [41]  and  [44].  It  is 
obvious  that  Equation  [41]  will  usually  apply  to  bulkheads  and  very  heavy  frames,  while 
Equation  [44]  will  apply  to  moderately  heavy  frames.  The  transition  point  will  vary  depending 
on  the  geometry  of  the  cylinders  and  the  lighter-frames.  Having  f|,  is  obtained  from 

Equations  [33]  and  [39]  as; 


For  the  rigid  bulkhead* 


•  is  allowed  to  go  to  aero  in  Equations  [41]  and  [45].  Some- 


oi^A p  4  5A)q 

times,  where  a  nearly  rigid  bulkhead  is  present,  the  actual  radial  deflection  of  the  bulkhead 
may  be  known  from  previous  experimental  data  more  accurately  than  an  effective  area  can  be 
computed.  Then  it  is  simply  a  matter  of  replacing 


.d[2«»/32pj(ij)  +  P4(l,)] 


by  yjdj)  in  Equation  [39].  Then  an  equation  analogous  to  Equation  [41]  results; 


} 


<f>9  «i)- 


^8(1  -  v*)  /  i\  r  /r  -  r  /  /I  X 

2(2-v)  V  r)  L“\r  + L  "  «  ‘  ‘  y.  J 


*-ff  h  (h)  -  0 

and  Equation  [45]  is  rsplaced  by 

^Ap 


i4«] 


- 1 


[47] 
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INTERSECTIONS 


At  an  axi symmetric  intersection  of  a  cylinder  with  another  shell  such  as  a  transition 
cone  or  an  ellipsoidal  or  hemispherical  eap,  the  simplest  method  of  design  is  to  estimate  the 
deflection  at  the  intersection  by  some  approximate  analysis  or  to  compare  the  proposed  struc* 
ture  with  available  experimental  results.  Then  Equations  [46]  and  [471  can  be  used,  or,  tor 
relatively  large  positive  (inward)  deflections.  Equation  [48]  should  be  replaced  by  a  modi  Red 
Equation  [44]. 


^4 


H  (h) 

if  I 


yi 


1 


[48] 


Of  course,  if  there  is  sufficient  knowledge  of  the  contiguous  shell  near  the  intersection 
at  high  pressure,  a  more  detailed  analysis  can  be  made  from  conditions  of  continuity  and 
equilibrium  at  the  intersection.  However,  this  will  usually  involve  more  work  than  is  practical. 


f 


-I 


■I 
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APPENDIX  B 


NUMERICAL  EXAMPLE 

A  Bumerioal  example  is  provided  to  illvstrate  the  prooedare  Car  desigBing  an  end  bay 
and  frame.  The  pertinent  scantlings  and  material  properties  are 

R  -  8.4180  in.  .  Rp  -  8.8464  in. 

h  m  0.0858  in.  v  -  0.80 

I  •  1.8217  in.  £  -  8  X  lO’  psi 

b  -  0.0448  in.  ffy  -  86,400  psi 

Aff  m  0.0516  in. 

The  cylinder  is  closed  by  a  rigid  bulkhead  ■  0],  and  is  estimated  to  be  1000  psi.  The 

computation  is  facilitated  by  use  of  an  instruction  sheet.  The  ftrst  column  of  the  instruction 
sheet  contains  the  item  number  in  parentheses.  The  second  column  indicates  the  operation, 
the  result  of  which  is  to  be  recorded  in  die  Uiird  column.  Numbers  in  parentheses  in  the  second 
column  are  item  numbers  and  indicate  that  the  operation  is  to  be  performed  on  the  quantity 
recorded  in  that  item,  e.g.,  (2)  (84)  means  multiply  the  quantity  in  Item  2  by  the  quantity  in 
Item  34,  whereas  2(34)  means  multiply  the  number  2  by  the  quantity  in  Item  84.  The  first  quan¬ 
tity  in  Item  84  is  always  -  0.1  depending  upon  the  sign  of  ^3(1),  Item  88.  The  subsequent 
values  are  entered  from  the  previous  interpolation  in  Item  121.  The  last  value  in  Item  86  is 
the  correct  value  of  1|,  and  Item  126  oont  .Ins  the  ratio  Ap  .  Thus  Lj,  the  distance 

from  the  bulkhead  to  the  center  of  the  first  frame,  is  obtained  from  by  adding  5/2,  and  Ap  ^ 
is  obtained  by  multiplying  Ap  by  Item  126.  In  the  example  which  follows  the  interpolation 
was  carried  out  three  times  in  order  to  demonstrate  that  the  first  interpolation  is  sufficiently 
accurate  for  design  purposes. 
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End  Bay  Desifiii  Calculation  Varying  and 


IlM 

Opcntion- 

R«mII 

1  ItM 

Operation 

Result 

f 

R 

1,4180 

(44) 

(lS)(43)/(7) 

0.265078 

(2) 

h 

0.0858 

(45) 

l-(44) 

0.734922 

(3) 

1 

L3217 

(46) 

l+(44) 

1.265078 

(4) 

b 

i? 

0.0443 

(47) 

^(457 

0.857276 

mi 

(S) 

Ag. _ 

'  Rm 

0.0491 

(48) 

V(46) 

L1247S7 

(6) 

0.3 

(49) 

(47) (48) 

0.964227 

(7) 

E 

3x  102 

(SO) 

(3)  (42) 

1.99906 

(8) 

-r 

66,400 

(51) 

(47)(50) 

1.7137 

(9) 

P. 

1,000 

(52) 

(48)  (50) 

2.2485 

[(10) 

VlUl) 

0 

(53) 

cosh  (51) 

2.86483 

(11) 

(6)* 

0.09 

(54) 

tinh  (51) 

2.68463 

(12) 

1-(S)*(11) 

0.79 

(55) 

cos  (52) 

-0.62701 

(13) 

1-2(6) 

0.4 

(56) 

sin  (52) 

0.77901 

(14) 

3-3(11) 

2.73 

(57) 

(48) (54) 

3.01969 

(IS) 

»'(i4r 

}.6S227 

(58) 

(47) (56) 

0.66785 

(IS) 

2-(6) 

1.7 

(59) 

(57)+(58) 

3.68754 

(17) 

(7)(10)/(1) 

0 

(60) 

(58)-(57) 

-2.35184 

(IS) 

m 

0 

(61) 

(53)-(55) 

3.49191 

(19) 

(13)(17)/|12) 

0 

(62) 

(49)  (61) 

3.36699 

(20) 

(8)2 

44J)896k  Vfl 

(63) 

(2) (4) 

0.0038 

(21) 

(12) (20) 

34.8a08x  to* 

(64) 

(5)+(63) 

0.0529 

(22) 

4(21)- 3  (18) 

139.3232xl0* 

(65) 

2{2)/(42) 

0.11346 

(23) 

V(ni 

11.803SX  19* 

(66) 

(S9)(64)+(6S)(62) 

0.57709 

(24) 

(23)/(12) 

14.9411x104 

(67) 

(5)/(66) 

0.085082 

(25) 

(24)+ (19) 

14.9411x104 

(68) 

(4)(42) 

0.067003 

(26) 

(l)/(2) 

98.112 

(69) 

(53)(56)/(48) 

1.98434 

(27) 

(9) (26) 

9.8112X  104 

(70) 

(54)(55)/(47) 

- 1.96349 

(28) 

(25)/(27) 

1.52286 

(71) 

(53)(55) 

-  L73623 

(79) 

(28)- 1 

0.52286 

(72) 

(54) (56) 

2.09153 

(30) 

(28)+  1 

2.52266 

(73) 

(68)(72)/(49) 

0.14534 

(3?) 

(29) (30)/ (28) 

0.86620 

(74) 

(44)(72)/(47) 

0.64672 

(32) 

(19)/(24) 

0 

(75) 

(71)-(74) 

-2.44295 

(33) 

(29)(32)/(30) 

0 

(76) 

(68) (75) 

-0.16368 

(34) 

l-(33) 

1 

(77) 

(69)- (70)+ (73) 

4.09317 

(35) 

(31)(34) 

0.86620 

(78) 

(69)+(70)+(76) 

-0.14283 

(36) 

(15)/(16) 

0.97192 

(79) 

(59) (67) 

0.31374 

(37) 

(35)(36)/2 

0.42094 

(80) 

(60) (67) 

-0.20010 

(38) 

(17)/ (27) 

0 

(81) 

l-(39)-(79) 

0.68626 

(39) 

2(38)/(16) 

0 

(82) 

(78)(81)/(77) 

-0.02395 

(40) 

(1)(2) 

0.722264 

(83) 

(37)+(8fl)+(82) 

0.19689 

(41) 

(15)/(40) 

2.28763 

(42) 

V?4i) 

1.51249 

(43) 

(26)(27)/2 

4.81290x104  1 

•DM* 

□ 

OpMlItM 

RomK* 

(*4) 

(121) 

-0.1 

-040974 

-047035 

-047033 

(») 

1-(I4) 

LI 

14CS74 

147035 

147033 

(M) 

(3)  (15) 

- 

1.4139 

14147 

1.4147 

'(•7) 

(51)  (15) 

1J»1 

LI332 

14343 

1.1342 

IB 

(51) (14) 

-0.1714 

-0.1195 

-0.1206 

-0.1205 

IS 

(52) (IS) 

2^734 

2.4053 

2.4007 

2.4060 

13 

(52)  (M) 

-0J24S 

-0.1501 

-0.1512 

-O.lSIl 

(SI) 

e«sh(l7) 

34114 

3.2009 

34102 

3.2099 

(12) 

sink  (17) 

34171 

3.0470 

34505 

34502 

(S3) 

cosh  (n) 

141472 

1.00715 

140721 

140727 

IB 

tinh  (II) 

-0.17224 

-0.11971 

-0.12019 

-0.12079 

IS 

cot  (IS) 

-0.71494 

-0.74096 

-0.74190 

is 

till  (SS) 

041957 

0.67155 

0.67051 

(ST) 

coo  (SO) 

047412 

0.91773 

041751 

0.91753 

(SI) 

tio  (SO) 

-0.22301 

-0.15616 

-0.15754 

-0.15744 

(SS) 

(S1)(II)+(I3)(S«) 

-0.12272 

0.17550 

0.10906 

0.17010 

(100) 

(S2)(I7)+(S4)(S5) 

3.27171 

3.09137 

3.10209 

3.10177 

(101) 

(Sl)(St)/(4S) 

145002 

141471 

1.91372 

1.91377 

(102) 

(S2)(S5)/(47) 

-244010 

-2.03351 

-2.03995 

-2.63944 

(103) 

(S1)(S5) 

-2.04471 

-2.37011 

-2.31105 

-241120 

(104) 

(S2)(SI) 

1.SS351 

2.04621 

2445391 

244543 

(lOS) 

(ll)(104)/(4S) 

0.UK3 

0.14219 

0.14213 

0.14213 

xm 

(44)  (104)/ (47) 

0.01042 

0.63271 

0.03245 

- 

(107) 

mum 

-3.20120 

-3.00119 

-341410 

- 

(101) 

(ei)(i07) 

-0.21151 

-0.20160 

-0.20195 

- 

(lOS) 

(47)  (98) 

-0.10520 

0.15050 

0.14545 

0.14512 

(110) 

(41) (100) 

3.07996 

3.41491 

341910 

3.41174 

(111) 

(losvf(no) 

347470 

3.63541 

3.63455 

343456 

(112) 

(109)-(110) 

-3.71510 

-343441 

-3.34365 

- 

(113) 

(101)-(102)4^(105) 

4.94065 

4.69046 

4.09510 

4.69534 

(114) 

(101)4-(102)-»(10S) 

-1.30K9 

-0.92047 

-0.92111 

- 

(115) 

(IT)  (111) 

0.30415 

0.30931 

040923 

0.30924 

(111) 

(17)012) 

-0.21370 

-0.21446 

■- 

(117) 

U(3S)-(U5) 

0.6SO0S 

0.69077 

0.6H76 

(111) 

(I14)(117)/(ll3) 

Hup 

-0.13554 

-0.13054 

- 

(US) 

(37)+(U6)+(lli) 

-041541 

0.00170 

-040006* 

- 

(120) 

(S3)-01S) 

041230 

0.19519 

0.19695 

- 

(121) 

(I4)(83)/(I20) 

-0.00974 

-0.07035 

- 

(122) 

1-(7S) 

- 

0.61626 

0.61626 

0.61626 

(123) 

(122) (113) 

- 

341119 

342254 

3.22222 

(124) 

(l5)/2 

- 

0.05673 

0.05673 

045673 

(125) 

(117)(124)/(123) 

- 

0.01216 

041216 

(121) 

l4-(125)/5 

m 

1441 

1.241 
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prrcpHure  h»^  bpcn  apvcloppd  for  compuUlinp  axisymmotnc  procedure  has  been  developed  for  coitipuUling  axi.symmetfic 

strp.^fiej  rear  these  disrmtinoitir.s  in  a  cylinder  in  which  the  stresses  near  these  discontinuities  in  a  cylinder  in  which  the 

thickness  of  the  raterial  does  not  vary  between  a  pair  of  stiff-  thickness  of  the  material  does  not  vary  between  a  pair  of  sUff- 

errrs  hut  ray  change  from  one  aide  of  a  stiffener  to  the  other.  pn,.rs  but  may  change  from  one  side  of  a  stiffener  to  the  other. 

The  results  of  this  analysis  indicate  that  stresses  higher  than  The  results  of  this  analysis  indicate  that  stresses  higher  tfian 


wdicwsl  by  t  sUndirf  compuinion  wi«  ajuMly  exist  new  these  discoelinuiUes.  A  those  t*edicted  by  s  sUndsri  compotstion  will  ususlly  exist  nesr  these  dis  .lU-  t i,  ..  \ 

relfcod  of  redacing  Uiis  effect  by  modifying  the  eeor.  try  near  these  poinU  is  also  presented.  method  of  redacing  this  effect  by  modifying  the  geometry  near  these  points  also  presented. 


